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INTRODUCTION 

Gas  absorption  is  one  of  the  numerous  chemical  engineering  operations 
which  involve  transfer  of  rsuterial  from  one  ohase  to  another.  Up  to  the 
present  time  nearly  all  industrial  absorption  equipment  is  operated 
under  the  cotintercurrent  flow  condition.  However,  evidence  indicates  that 
in  certain  cases  this  type  of  flow  is  not  necessarily  the  best  (17). 

The  main  disadvantage  of  co-current  flow  is  that  generally  a  smaller 
driving  force  is  established  in  the  absorption  equipment  in  comparison 
with  the  countercurrent  flow.  Such  disadvantage  will  disappear  when  there 
is  no  appreciable  vapor  pressure  of  the  transferring  component  in  the 
liquid  phase.  For  instance,  in  the  absorption  of  CO,  into  caustic  solutions, 
it  has  been  assumed  by  various  investigators  (17),  that  there  is  no  COg 
pressure  over  the  solution.  Thus,  the  driving  forces  in  both  countercurrent 
and  co-current  flows  will  be  the  same,  as  shown  by  the  following  graph. 
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On  the  other  hand,  operation  under  co-current  flow  has  some  advantages 
over  that  under  countercurrent  flow.  The  pressure  drop  is  smaller  in  the 


former  case,  and  correspondingly  the  power  requirement  will  be  lower. 
Besides,  there  is  no  limitation  of  high  rates  of  flow  due  to  flooding  or 
loading. 

Various  apparatus  have  been  proposed  for  laboratory  scale  absorption 
studies.  Two  of  the  frequently  used  apparatus  are  the  disc  column,  which 
was  introduced  by  Stephens  and  Morris  (19),  and  the  short  wetted-wall 
column,  which  has  been  widely  used  during  the  last  few  years  (4.)  (22). 
In  both  columns  the  surface  areas  are  measurable,  but  the  disc  column 
provides  a  condition  more  like  that  in  the  packed  absorber. 

This  investigation  was  to  study  the  co-current  flow  absorption.  The 
experiments  were  performed  in  both  a  disc  and  a  packed  column  under  co- 
current  and  counter-current  flows  for  the  sake  of  comparison.  The  C02-H20 
system  was  chosen  for  this  study,  not  only  for  the  well-known  fact  (16) 
that  the  liquid-film  is  the  only  controlling  factor,  but  also  for  the 
extensive  results  of  other  investigators  available  for  comparison  (6). 

THEORI 

In  developing  absorption  theories,  several  models  have  been  proposed 
to  describe  the  absorption  mechanism. 

Hatta's  Film  Theory  (15) 

The  film  concept  pictures  a  stagnant  fluid  film  at  the  interface, 
through  which  film  the  substance  to  be  absorbed  is  transferred  by 
stationary  diffusion.  The  rest  of  the  liquid  is  considered  to  be 
completely  homogenized. 


Surface  Renewal  Theory  (3) 

This  is  the  modified  form  of  the  Higbie  Penetration  Theory  due  to 
Danckwerts.  In  this  theory  an  element  of  the  liquid  present  at  the  inter- 
face is  changed  by  a  transient  diffusion  process.  After  some  time  the 
element  is  replaced  by  another.  The  char.ce  of  the  element  being  replaced 
within  a  given  time  is  assumed  to  be  a  statistical  distribution  and 
independent  of  its  age. 

In  both  concepts  use  is  made  of  a  quantity  which  can  not  be  directly 
measured.  In  the  film  theory  this  is  the  effective  film  thickness  x., 
while  in  Danckwerts'  theory  it  is  the  mean  rate  of  production  of  fresh 
surface  per  unit  surface  F.  According  to  these  theories,  the  liquid  film 
coefficient  can  be  expressed  as: 

Tt^=  D/xj  by  Hatta*  s  film  theory,  and 
Tl^—JISF    by  surface  renewal  theory, 
where  kL  is  the  liquid  film  coefficient  of  physical  mass  transfer,  and  0 
is  the  diffusivity. 

Boundary  Layer  Theory  (14) 

In  this  theory  the  diffusion  boundary  layer  is  considered.  The 
thickness  of  the  layer  is  the  distance  measured  normally  from  the  interface 
in  which  the  concentration  changes  from  the  interfacial  value  to  the  stream 
value.  This  theory  differs  from  the  film  theory  on  the  point  that  the 
boundary  layer  is  moving  while  the  film  is  assumed  stagnant.  This  theory 
was  developed  only  recently.  Owing  to  mathematical  difficulties,  only  some 
simple  cases  have  been  treated.  However,  the  theory  provided  a  more 
realistic  physical  picture  than  any  other  theory. 


The  applicability  of  the  boundary  layer  theory  will  likely  be  reduced 
considerably  by  the  instability  of  a  fluid-fluid  interface,  but  the  same 
is  true  of  the  stagnant  film  and  Surface  Renewal  theories.  Some  Important 
papers  on  this  field  are  listed  in  the  references  (12)  (H). 

EXPERIMENTAL 
Scope 

Two  types  of  columns,  disc  column  and  packed  column,  were  used  in  the 
present  investigation.  In  the  disc  column,  the  range  of  liquid  flow  was 
determined  by  the  rates  at  which  the  column  became  unstable;  that  is,  from 
the  lowest  rate  which  maintained  nearly  perfect  wetting  (cf .  section  under 
the  heading  'De-wetting*  on  p.  14-)  of  the  disc  surface,  10.0  lb./hr. ,  to 
the  highest  rate  which  kept  the  water  from  spraying-out,  33.0  lb./hr.  In 
the  packed  column,  the  liquid  rate  covered  a  range  of  13.0-52.6  lb./hr. 
The  gas  rate  varied  from  1.55-10.6  cu.ft./hr.  (the  corresponding  Reynolds 
numbers  based  on  effective  column  diameter  were  606  and  4150  respectively). 
The  highest  liquid  temperature  was  31.6°C,  and  the  lowest  was  21.0°C. 
The  column  pressures  fluctuated  between  722  -  755  mm.  Hg.  The  experimental 
quantity  determined  was  the  liquid  film  coefficient  at  25°C. 

Equipment 

Columns,  The  disc  column  consisted  of  35  eeranic  discs,  enclosed  in  a 
pyrex  glass  tube  of  1-1/8  inch  inside  diameter.  The  discs  were  threaded 
edgewise  on  a  vertical  fiberglas  cord  in  such  a  way  that  the  successive 
ones  were  maintained  at  right  angles  by  means  of  Duco  cement. 


The  general  arrangement  is  shown  on  Plate  I,  The  water  was  introduced 
at  the  top  through  a  jet,  and  removed  by  a  central  tube  and  small  funnel 
under  the  lowest  disc.  The  liquid  feed  jet  was  placed  5  cm.  above  the 
uppermost  disc,  as  recommended  by  Stephens  and  Morris  (19). 

The  packed  column  was  constructed  with  a  2  inch  inside  diameter  pyrex 
glass  pipe,  packed  with  8  mm.  glass  rasching  rings.  The  bed  was  5-2/4" 
in  height  and  supported  by  a  perforated  plate.  Water  was  distributed  over 
the  packings  by  a  Tygon  sprayer.  Both  ends  of  column  were  connected  to 
glass  tees,  which  formed  the  gas  calming  sections.  Other  constants  for 
both  columns  are  listed  in  Appendix  I. 

Accessories.  A  1/8  hp.  centrifugal  pump  was  used  to  feed  water  (distilled 
water)  from  a  5  gal.  carboy  to  tlie  absorption  columns.  In  the  case  of  the 
disc  column,  a  constant  head  tank  was  used. 

A  1/4"  needle  valve  was  placed  before  a  f lowrator  to  regulate  the 
liquid  flow. 

Another  1/4"  needle  valve  was  placed  before  a  rotameter  to  regulate 
the  gas  flow  from  a  COg  cylinder. 

Gas  leaving  the  rotameter  passed  through  three  saturation  bottles. 
The  difference  between  the  inlet  and  outlet  gas  temperature  was  kept 
within  1°F  to  eliminate  any  effect  due  to  vaporization  of  water.  The 
guage  pressure  of  the  gas  flow  was  measured  by  a  manometer. 

A  sampling  reservoir  was  used  to  stabilize  the  outlet  liquid  flow 
during  sampling. 

Themometers  were  provided  at  the  inlets  and  outlets  of  both  gas  and 
liquid  lines. 
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Procedure 


1.  Discharge  remaining  liquid  from  the  column  and  sampling  reservoir. 

2.  Set  gas  rate  G  and  liquid  rate  L. 

3.  Adjust  the  height  of  the  sampling  reservoir  to  ensure  liquid  seal 
with  minimum  liquid  level  in  the  liquid  collecting  tube. 

A*       Take  the  first  sample  at  10  min.  for  the  disc  column,  or  at  5  min. 
for  the  packed  column  after  the  start  of  the  experiments. 

5.  Sampling  tubs  was  dipped  into  the  trapping  solution  (20  ml.  of  0.1H 
BaCl2  and  40  ml.  of  0.05N  NaOH)  in  a  100  ml.  measuring  flask.  The 
flask  was  plugged  with  rubber  stopper  before  filtration. 

6.  Drain  off  the  liquid  remaining  in  the  sampling  tube. 

7.  Take  the  second  sample  at  16  min.  for  the  disc  column,  or  at  8-10 
■in.  for  the  packed  column. 

8.  Time  of  sampling  was  about  4.5". 

Chemical  Analysis 

Standard  solutions  t,  BC1  wa3  standardized  against  NaoG0,.  The  effective 
concentration  of  NaOH  (OJTion  onl7,  C03"not  considered)  was  standardized 
against  the  standard  HC1  solution,  after  excess  BaCl2  was  added. 
Inlet  H20.  The  concentration  of  003  in  the  inlet  distilled  water  was 
analysed  every  six  runs.  No  appreciable  content  of  CO2  in  H20  was  found. 
Purity  of  CO?.  The  purity  of  cylinder  C02  was  analyzed  by  volumetric 
method.  The  result  showed  that  it  contained  C02  wore  than  99.2  percent. 
In  calculation,  100  percent  was  assumed. 

Analysis  of  samples.  The  absorbed  003  was  precipitated  as  BaCOj  In  the 
trapping  solution.  After  filtration  the  excess  NaOH  was  determined  by 


back  titration  against  HC1  solution  with  phenolphthalein  as  indicator. 

A  Magnetic  stirrer  was  used  in  titration. 

Filtration  had  no  appreciable  effect  on  the  analytical  results,  as 
shown  in  Appendix  II. 

Preliminary  Experiraents 

A  suspicion  of  errors  introduced  by  the  filtration  of  the  sampling 
solution  in  the  course  of  analysis  led  to  the  conduction  of  an  auxiliary 
experiment.  Its  results  are  summarized  in  Appendix  II.  It  was  shown  in 
this  auxiliary  experiment,  that  the  results  of  analysis  were  practically 
unaffected  by  filtration,  the  variation  of  the  height  of  funnel  nozzle 
above  filtrate  surface,  or  the  exposure  time  when  it  was  less  than  30 
minutes. 

The  same  experiment  showed  that  the  variation  of  sample  obtained  in  the 
disc  column  was  within  experimental  error  for  samples  taken  at  longer  than 
8  minutes  after  start  of  the  experiment,  and  another  experiment  showed 
that  in  the  packed  column  it  was  A  minutes  after  start  (Appendix  II ). 

Design  of  Experiment 

In  the  design  of  experiment  letters  P,Q,  and  R  denote  the  co-current 
flow,  and  C,D,  and  E  the  counter-current  flow  operation. 

A  2x3x3  factorial  experiment,  PC  (2  types  of  flow,  3  levels  of  gas 
rate  and  3  levels  of  liquid  rate)  was  conducted  for  the  disc  column  study. 
This  experiment  contained  2  replicates  with  total  36  runs  and  72  observations. 

Another  2x5  factorial  experiment,  QD  (2  types  of  flow  and  5  levels  of 
liquid  rate)  with  various  repititions  was  conducted  for  the  disc  column  af 
a  fixed  gas  rate. 
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A  2x3x4.  factorial  experiment,  HE  (2  types  of  flow,  3   levels  of  gas 
rate  and  U  levels  of  liquid  rote)  with  2  replicates  was  conducted  for  the 
packed  column  study. 

Randomized  complete  block  designs  were  used  for  the  experiments  PC  and 
RE,  and  completely  randomized  design  was  used  for  the  experiment  QD. 

RESULT  AND  DISCUSSION  OP  EXPERIMENT 

Calculation  of  Liquid-Film  Coefficients 

For  calculating  liquid-film  coefficients  from  the  experimental  data, 

the  following  equations  were  used ; 

for  the  disc  column  kr  ~   .  R^*  (l) 

j.»m* 

for  the  packed  coluan  ,  waj  ,  . 

kL*  --  ,  C  (2) 

<*  C>l.m. 

where 

kL=  liquid-film  coefficient,  lb./(hr.-sq.ft.) (lb./cu.ft.) 

ki£=liquid-fiLn  coefficient  on  a  volume  basis,  lb./(hr.-cu.  ft.)  (lb./cu.ft.) 

H  —  rate  of  absorption,  lb./hr. 

V  -   volume  of  the  packed  bed,  cu.ft. 

A  =  dry  surface  area  of  the  discs,  3q.ft. 

^C)l.«u=  logarithmic  mean  of  (Ce-C)  at  inlet  and  outlet,  lb./cu.ft. 

Ce  =  liquid  concentration  in  equilibrium  with  the  gas  phase,  lb./cu.ft. 

C  =  liquid  concentration,  lb./cu.ft. 
The  values  of  Ce  were  calculated  from  the  Henry's  law,  Ce  =p/H,  where  p  ig 
the  partial  pressure  of  COg,  and  H  the  Henry  constant  (13).  The  reason 
for  replacing  the  interfacial  concentration  with  the  equilibrium  concentration 
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Cft  is  based  on  the  fact  that  the  CO.- -HO  absorption  is  controlled  by  the 
liquid  film,  as  verified  by  Sherwood  and  Kolloway(l6). 

All  mss  transfer  coefficients  so  obtained  were  corrected  to  25°C 
according  to  the  following  equation,  (16): 

kL5  -  kL.e°-0a3(25~t)     (3) 
Summary  of  Data 

Data  are  suoitarized  on  Table  3  to  7  in  Appendix  III.  Table  3  lists 
observed  data  for  PC  series  experiments.  Table  A  gives  values  of  liquid- 
film  coefficients  calculated  from  the  data  in  Table  3.  Table  5  contains 
both  observed  data  and  calculated  liquid-film  coefficients  for  QD  series 
experiments.  Table  6  contains  observed  data  for  RE  series  experiments 
and  Table  7  lists  the  values  of  liquid-film  coefficients  calculated  from 
the  data  in  Table  6. 

Analysis  of  Data 

Data  obtained  from  experiments  PC  and  RE  were  analyzed  statistically. 
The  results  are  given  in  Table  1  and  2  below.  The  detail  can  be  found 
in  Appendix  IV  (p.  £i  ).  It  is  Irnportant  to  notice  that  the  ordinarily 
assuwd  additive  model  In  statistical  analysis  is  not  applicable  to  the 
kL  value,  since  multiplication  is  involved  in  the  evaluation  of  kT. 
Such  operation  will  lead  to  serious  error  in  standard  deviation,  and  thus 
transformation  is  necessary  before  analysis.  The  logarithmic  transforma- 
tion was  carried  out  for  this  purpose. 
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Table  1.  Analysis  of 

variance]  FO 

'-3S. 

Source  of  ..yaci&'fciQfl  ., 

d.f. 

S.3- 

Bad 

Replicate  (R) 

1 

0.19995 

0.19995 

Treatment  (T) 

17 

2.04449 

0.12030 

I'low  tyne  (F) 

1 

0.02467 

0.02467 

liquid  rate  (L) 

2 

1.89697 

0.94850 

Gas  rate  (G) 

2 

0.02553 

0.01276 

FxL 

2 

0.00241 

0.00121 

FxG 

2 

0.01207 

0.00604 

I«fl 

4 

0.00962 

0.00240 

FxLxG 

4 

0.07322 

0.01330 

Error 

53 

0.24888 

0.00470 

RxT 

17 

0.13487 

0.00793 

Obs'n  :  run 

36 

0.1L401 

0.00316 

Total 

71 

2.49332 

Table  2.  Analysis  of 

variance;  RE 

1 pries. 

Source  of  Variation 

m.s. 

Replicate  (R) 

1 

0.00206 

0.00206 

Treatment  (T) 

23 

3.65824 

0.15905 

Flow  type  (F) 

1 

%  00270 

0.00270 

Liquid  rate  (L) 

3 

3.62859 

1.20963 

Gas  rate  (G) 

2 

0.00079 

0.00040 

FxL 

3 

O.OO403 

0.00134 

FxG 

2 

0.0^552 

0.00276 

LxG 

6 

0.00472 

0.00079 

FxLxG 

6 

0.01189 

0.00198 

Error 

69 

0.06956 

0.00101 

RxT 

23 

9.0422? 

0.00183 

Obs'n  :  run 

46 

0.02729 

0.000593 

Total 

93 

3.72986 

Since  the  m.s.  of  the  main  affects  (F,L,  and  G)  estimate  error  terns  of 

very  complicated  forms  (5), 

it  is  not  suitable  to  use  Snodoeor's  F  to  test 

the  significance  of  those  effects.  However,  comparing  the  main  effects  with 

interaction  terms  (FxL  etc. 

),  it  is  reasonable  to  believe 

that  there  is  no 

difference  between  flow  types,  and  only 

liquid  rate  has  any  effect  on  the 

liquid-film  coefficient  within  the  range 

of  study. 
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The  significance  of  variation  due  to  replicate  in  PC  series  (Table  l) 
indicates  that  soiae  unnoticed  error  might  be  introduced  in  either  of  the 
two  replicates.  A  discussion  on  this  case  is  given  in  Appendix  V  under 
the  heading,  Selection  of  Data  (p. C^ )• 

The  standard  deviation  is  13.8  percent  for  PC  series,  A  percent  for 
QD  series,  and  5.8  percent  for  RE  series.  The  high  deviation  in  PC 
series  must  result  from  a  few  scattered  data. 

Experimental  Error 

Error  Due  to  Liquid  Sampling.  In  order  to  get  better  reading  of  sample 
volume,  the  sampling  nossle  was  immersed  le3S  than  1-1/2  inches  below  the 
surface  of  liquid  in  the  measuring  flask.  But  this  could  not  ensure  no 
loss  of  gas  from  samples.  The  magnitude  of  such  error  is  not  easy  to 
estimate. 

Error  Due  to  Method  of  Analysis.  The  method  of  Emmert  and  Pi^ford  (A)  was 
used.  An  average  error  of  2  pereent  was  supposed  to  be  involved.  As 
MBticned  by  Taylor  and  Roberts  (20),  this  method  was  batter  than  the 
barium  hydroxide  method,  which  generally  results  3  percent  in  error,  and 
the  method  of  Uaanerton  and  Garner,  which  gives  results  systematically 
hifh  by  about  10  percent. 

Error  Diis   to  Liquid  Temperature.  Thermometers  with  scale  graded  to  0.1°C 
were  used  on  liquid  lines  for  QD  and  HE  series,  and  with  scale  graded  to 
0.5*^  were  used  for  series  PC.  A  misreading  in  0.1°C  would  lead  to  an 
error  in  equilibrium  constant  of  about  0.3  oercent,  and  this  would  further 
be  enlarged  by  the  k^  computation  formula  to  about  0.9  percent,  as  will  be 
whown  in  a  later  section.  In  RE  series  experiments,  the  pump  caused  the 
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inlet  water  tempera ture  to  continuously  increase  at  a  low  flow  rate.  This 
complicated  the  absorption  process  by  sensible  heat  transfer.  This  un- 
steadiness of  the  absorption  condition  would  give  rise  to  some  deviation, 
and  this  was  counted  as  experimental  error.  The  total  error  contributed 
by  liquid  temperature  deviation  to  the  value  of  kL  is  thus  considered  as 
1.5  percent. 

De-wetting  of  the  Liquid  on  the  Disc.  The  de-wetting  phenomenon  has  been 
reported  by  a  number  of  investigators,  and  in  some  cases  the  de-wetted 
areas  have  been  quite  extensive,  even  at  liquid  flow  rates  up  to  200 
Ib./hr.-ft.  (20).  In  this  experiment,  de-wetting  was  found  at  a  liquid 
rate  as  high  as  173  lb./hr.ft.,  e.g.  run  Q-9-4,  and  complete  wetting  was 
found  at  the  liquid  rate  as  low  as  103  lb./hr.ft.,  e.g.  run  D-5-1. 
Generally,  de-wetting  rarely  occurred  at  the  rate  higher  than  120  lb./hr.ft. , 
and  complete  wetting  was  hard  to  find  at  the  rate  below  129  lb./hr.ft.  The 
largest  de-wetted  area  observed  visually  was  about  14.  percent  of  the  total 
area.  Loss  of  absorption  surface  will  cause  low  absorption  coefficients, 
while  the  increased  flow  rate  in  the  wetted  areas  will  tend  to  offset  this. 
It  is  possible  that  partial  de-wetting  might  give  rise  to  either  high  or 
low  results. 

It  was  found  that,  at  high  flow  rate,  liquid  would  sometimes  drop  from 
one  disc  to  the  next  paralleled  to  it,  without  touching  the  neighboring  one, 
which  was  at  tb*  right  angle  to  it.  This  phenomenon  would  also  decrease  the 
contact  surface  area  of  liquid  and  gas,  and  resulted  a  relatively  low  co- 
efficient. 

Errors  Amplified  bv  Transfer  Enuation.  Errors  will  be  amplified  2  or  3 
times  by  the  transfer  equation.  This  will  be  shown  by  the  following 
illustration: 
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Illustration:    Run  POO-1  had  the  following  observed  data: 

p.p.  of  CO  ~  column  pressure  -  p.p.  of  H20  —  730-24.  -  706  mm. 
Liquid  temperature  both  at  inlet  and  outlet  ^24.5°C. 
The  Henry's  Law  constant  at  24.5°C  is  8.1QX163  mm.Hg/lb.COycu.ft.   The 
inlet  concentration  of  C02  in  H20  is  zero,  and  that  at  the  outlet  is 
0.0470  lb./cu.ft.  kL  was  computed  accouding  to  eq.  (l),  where  •  «=  CxL/62.4. 
(L  is  the  liquid  mass  velocity  in  unit  of  lb./hr.  per  area  of  contact 
surface  in  sq.ft.).  The  value  of  k  -  0.566  at  24.5°C  was  obtained. 
If  a  deviation  of  5  percent  less  than  the  present  value  of  C  occurs,  kT 
will  be  equal  to  0.517,  i.e.  it  causes  an  error  of  9.15  percent,  and  if 
a  deviation  of  0.5  percent  less  than  the  present  value  of  CQ  occurs,  kr 
will  become  0.574,  and  an  error  of  1.4  percent  results.  In  the  case  of 
different  Ce  at  inlet  and  outlet  the  deviation  will  be  even  larger. 

DISCUSSION  OF  RESULTS 

Effect  of  Flew  Type 

It  was  expected  that  co-current  flow  and  countercurrent  flow  would 
have  some  different  effects  on  the  mass  transfer  coefficient.  According 
to  the  stagnant  film  concept,  the  liquid  film  should  be  thinned  and  the 
holding  time  should  be  much  decreased  by  the  co-current  flow.  According 
to  the  boundary  layer  theory,  the  relative  velocity  should  play  an  important 
part  in  establishing  the  thickness  of  fluid  layer.  However,  this  inves- 
tigation showed  that  there  is  no  effect  of  using  co-current  flow  or 
countercurrent  flow  on  kL  within  the  range  of  investigation.  In  studying 
the  COg-ajO  absorption  system  in  wetted-wall  columns,  Collins  (2)  found 
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that  the  use  of  co-current  flow,  at  Reynolds  number  of  gas  higher  than 
14,800,  increased  the  transfer  coefficient  appreciably. 

Effect  of  Gas  Rate 

In  this  investigation,  the  gas  rate  range  was  too  .narrow  to  detect 
any  influence  on  the  liquid  film  coefficient.  Hildta  et  al.  (7)  found 
that  the  liquid-film  coefficient  of  CO^-RgO  system  in  a  wetted-wall  column 
was  affected  by  gas  rates  at  Reynolds  number  greater  than  6000  in  counter- 
current  flow,  when  the  liquid  rate  was  such  that  Re  -  300  and  also  Re-  600. 
This  seems  contradictory  to  what  might  be  expected  by  boundary  layer  theory 
for  the  simplest  case  (2-dimensional,  co-current  flow  with  horizontal 
interface),  in  which  the  relative  velocity  as  well  as  the  absolute  velocity 
of  gas  is  the  determining  factor  for  rate  effect  on  transfer  coefficient. 
It  could  be  explained  as  that  the  effect  of  gas  rate  was  due  to  ripple 
formation  rather  than  the  change  in  film  thickness. 

Effect  of  Liquid  Rate 

Since  there  was  no  interaction  between  "lew  type  and  liquid  rate,  as 
shwm  by  analysis  (Table  1  and  2),  a  single  correlation  shall  be  provided 
for  each  column.  Logarithmic  plots  of  kL  vs.  p  and  kLa  vs.  L  were  con- 
structed for  the  disc  column  and  packed  column  respectively.   pis  the 
wetting  rate  (equal  to  the  liquid  flow  rate  in  lb./hr.  divided  by  the  mean 
perimeter  for  liquid  flow  in  ft.)  and  L  is  the  liquid  flow  rate  In 
lb./hr. sq.ft.  These  lines  can  be  represented  in  the  following  form: 

or  „ 

kLa  -=  bLn 


where  b  and  n  are  the  constants  to  be  determine:'-  experimentally. 

For  the  disc  column  the  absorption  data  can  be  correlated  in  the 
following  equation: 

k^  0.0203  P  °'7A5        »*  20°C     U) 
The  equation  was  based  on  88  observations  (Appendix  V).  The  sample 
standard  deviation  from  the  equation  is  0.0436  in  logarithmic  scale,  or 
10.5  percent  of  the  value  of  kj^.  The  sample  standard  deviation  of  the 
slope  is  0.0298,  or  4-.0  percent. 

For  the  packed  column  the  absorption  data  can  be  correlated  by  the 
equation: 

kLa  -  0.655  L0,85    at  20°C    (5) 
The  equation  was  based  on  88  observations.  The  sample  standard  deviation 
from  regression  is  0.0284.  in  logarithmic  scale,  or  6.75  percent  of  the 
value  of  kja.  The  sample  standard  deviation  of  the  slope  (the  regression 
coefficient)  is  0.0127,  or  1.5  percent. 

Generalized  Correlation 

Sherwood  and  Holloway  (16),  investigated  desorption  of  oxygen  from 
water  and  they  proposed  the  following  generalized  correlation  for  kT : 

D      -*       X 

*L*_4(-ir)n  (-ft)*         (7) 

where>.is  the  viscosity  of  the  liquid,  (>  is  the  density  of  the  liquid, 
and  c*i  and  n  are  constants  of  a  particular  column.  When  the  general 
equation  is  applied  to  the  results  of  this  experiment,  the  following 
equations  are  obtained: 
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for  the  disc  column        J^  V.U^)0*745-^)^     (8) 
for  the  pecked  column        kLa   8^.0(Jp-)0'85  (JL)05  (9) 

Comparison  with  Results  of  Previous  Workers 

Disc  column  equation.  The  liquid-film  coefficients  obtained  are  lower 
than  all  the  published  results.  Stephens  and  Morris  (19)  have  mentioned 
that  the  absorption  coefficients  obtained  on  different  disc  columns  might 
vary  by  ±  10  percent.  However,  data  with  deviation  about  50  percent 
lower  than  that  given  by  Stephens  and  Morris  have  been  found  in  the 
literature  (6).  The  present  result,  though  much  lower,  gives  a  line 
aearly  parallel  to  Stephen  and  Morris' .  Their  data  was  represented  by 
the  equation: 

i^  =  0.048  r  °'7 

In  plotting  their  data  for  GO2  absorption  in  Doulton  disc  and 

pyrophyllite  disc  columns,  Taylor  and  Roberts  (20)  observed  the  existence 

of  a  distinct  change  of  slope  in  liquid  film  coefficient  versus  wetting 

rate  plot.  Their  results  for  both  columns  were  correlated  into  a  single 

set  of  equations,  via: 

For      H<  1551b.  Ar.ft.  k     =  0.124  V°'U 

L 

For   T>155  Ib./hr.ft.  kj  =  O.OO&r1*0 

Chu  (l),  using  the  same  column  as  that  for  the  presenc  study,  also  observed 
the  break  slopes  in  his  data.  He  obtained  the  following  set  of  equations: 
For   73  <P< 200  lb./hr.ft.       k.  =  O.^Son0*6^    (10) 
For   16  <p<  73  lb./hr.ft.      kL  m   0.123  T  °'37 
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The  difference  in  critical  flow  rates  at  which  the  break  occurred  has  been 

reported  by  Taylor  and  Roberts  (20)  after  the  study  of  six  different  types 

of  disc  columns.  However  no  such  break  was  observed  in  this  study,  ncr 

in  Hwu's  work  (9).  Hwu  constructed  the  present  column,  and  he  suggested 

the  equation  0  Q, 

^   -  0.0075  T  (H) 

for  GO  absorption  in  this  column.  The  liquid-film  coefficients  found  by 

Hwu  were  higher  than  those  of  Ghu  (l)  and  the  present  investigator.  All 

the  results  just  mentioned  and  some  others  are  plotted  on  Plate  III  for 

comparison. 

Prediction  of  Correlation  for  CI.7-H  0  System  by  the  Present  Result.  The 

absorption  of  chlorine  is  a  typical  liquid-film  controlled  system,  as  has 

been  shown  by  Sherwood  on  the  basis  of  CC*  and  oxygen  absorption  and  de— 

sorption  data.  It  has  also  been  recognized  that  in  the  tower  with  a  small 

diameter  the  variation  of  gas  rate  has  no  effect  on  the  transfer  coefficient, 

as  verified  by  Vivian  and  Whitney  (23).  Therefore,  the  result  obtained  from 

COg  absorption  study  is  expected  to  be  applicable  to  the  chlorine-water 

absorption  data. 

Using  the  general  equation  (6),  or  remembering  that  3c  varies  with 

0.5 
D  '  for  the  same  absorbent  at  the  same  liquid  rate,  we  can  derive  an 

equation  for  CLj-HjO  system  as: 

1^  =  0.0178  P0-7*5        (12) 

Similarly  we  have  the  corresponding  equations  derived  from  Chu's  equation 

and  Ewu's  equation.  These  are J 

kL  =  0.0338  P0,644        (13) 
kL=  0.00657/-' °'95         (U) 

respectively.  The  three  predicted  equations  are  represented  by  lines  II, 


I 

P 

a) 


P 


p 


5 


o 

CQ 

'O 

O 
I 


o 

u    • 
1  t? 

p  p 

©  c 


—».     o 


g 

1 

jq 

■ 

O 

Vl 

<w 

O 

o 

P 

P 

i 

i 

in 

CO 

1 

I    I 


9> 


O 


M 
H 

I 

a.. 
o 

SB 

o 

§ 


s 


cr 
a 
I 

• 

i 


o 

•H 

© 
O 
O 

A 

3 


I 

jr 

P 

P 


© 
p 


% 

a 

© 
p 

g 

p 
p 

I 


I  I 


© 

£ 
p 

© 

© 

I 


o^     o 


© 


3 
© 

•a 
© 
p 

CO 

«M 
O 

P 

i 

2 


© 
© 
I 

3     so 
o 


O 
P 

| 

© 


M 
M 


© 


o 

V( 

o 
p 

1 

© 


|      I      <      I 


O 


O 


o 


Ill,  and  IV  on  Plate  IV,  and  compared  with  the  line  directly  drawn  from 
the  experimental  data  by  Tien  (2l)  in  the  same  column.  The  corresponding 
equation  for  the  experimental  data  is: 

kL  -  0.0l63r0'8:i- 
The  agreement  of  equation  (12)  with  that  obtained  by  Tien  within  the  range 
of  experiment  (llO-^OO  lb./hr.i't.)  is  clearly  shown  by  Plate  IV.  The 
predicted  values  of  liquid  film  coefficient  by  use  of  Chu's  equation  (13) 
are  little  higher  than  the  experimental  values,  and  those  predicted  by 
Hwu*s  equation  were  even  higher.  The  derivation  of  the  equations,  and  their 
representative  points  are  given  in  Appendix  VII . 

Prediction  of  Gas-film  Coefficient  of  HHy-HgQ  system.  By  the  eosbination 
of  liquid-film  coefficient  data  and  overall  mass-transfer  coefficient,  K„, 
data,  we  can  calculate  the  gas-film  coefficient  from  the  following 
relationship: 

1    .  J  ,.  -   H 
*G       %       \ 

where  H  is  Kenry  constant. 

Hwu  has  determined  the  overall  mass  transfer  coefficient  of  NEU-H20 
system,  and  calculated  the  gas-film  coefficient  by  using  his  own  equation 
for  liquid-film  coefficient.  Since  equation  (12)  predicts  liquid  film 
coefficients  in  the  Cl2-H20  system  better  than  Hwu*s  equation,  an  attempt 
was  thus  made  to  use  equation  (^-)  of  the  present  investigation  together 
with  Hwu's  experimental  data  of  Kg  to  calculate  kg  for  M-J^O  system. 
The  detail  of  calculation  are  given  in  Appendix  VII.  The  resulting 
equation  Is: 

Xq  =3.90r°*3°  (15) 


EXPLANATION  OF  PLATE  IV 

Chlorine-water  absorption  correlation  (k^  vs.f1  ) 
in  the  disc  column 


k  liquid  film  coefficient,  lb./(hr.-sq.ft. ) (lb./cu.f t. )  at  20°C 

p  wetting  rate  (liquid  rate  oer  mean  wetting  peritieter  of  disc), 

lb./hr.ft. 

Curve  I  mm  experimental  result  by  Tien  (21) 

Curve  II  —  predicted  by  the  present  work  on  the  basis  of  GOj-HgO  system 

Curve  III  —  predicted  by  the  result  of  Chu*s  work  on  the  basis  of 
COjj-ELO  system 

Curve  IV  —  predicted  by  toe  result  of  Hva^s  work  on  the  basis  of 
COg-Hg0  systtB 
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or  in  the  general  form 

V  _  0.0326 rr-30(frO.33(A)-.56^)  (l6) 

The  equation  (15)  was  corrected  to  a  relative  velocity  of  8.4  ft./sec. 
for  the  convenience  of  comparison  with  published  data  (Hwu  corrected  his 
data  to  the  relative  velocity  5.84  ft./3ec).  Plate  V  shows  the  comparison 
of  the  results  from  various  sources.  It  may  be  noted  that,  over  the  range 
of   studied  by  Hwu  (  p=  155-395  Ib./hr.ft.),  equation  (15)  is  quite  con- 
sistent with  the  experiment  data  of  Taylor  and  Roberts.  Hwu's  equation 
corrected  to  8.4.  ft./sec.  gives  relatively  low  gas-film  coefficients.  The 
equations  of  the  curves  are: 
S.  andM.  (19)  ifc  ^ll.lp0,23 

T.  and  R.  (20)         fy  a    5.3  p0'25 
Hwu  (9)  kfc  =2.99  p0,26 

Packed  Column  Equation.  Koch  et  al.  (10)  studied  CO2  absorption  in  6-  and 
J0-ineh  towers  with  a  considerable  variety  of  packing  rings.  He  correlated 
all  his  data  by  the  equation. 

Kja  =  0.015  h°        ,  or 

HqL  =  1.05  L0#Q4 
where  the  result  of  the  present  investigation,  equation  (5)  and  the 
corresponding  equation: 

HoL  «  0.103  L0'15  (17) 

^HoL  "  ^Al®)  is  the  height  of  transfer  unit  ^ive  much  higher  liquid-film 
coefficients  at  the  low  liquid  rate  range  than  Kod^s. 

Since  the  diffusivity  of  CO.,  and  Oj  are  6.8-7.0  xlO**5  sq.ft./hr.  at 
20°C  (11)  (17),  equation  (5)  should  directly  be  applicable  to  O^j-H-O 
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system.  A  comparison  with  02-H20  absorption  data  by  Sherwood  and  Holloway 
(16)  is  shown  on  Plate  Tin.  The  present  results  give  higher  HoL  values. 

Ghu  (l)  worked  with  the  same  CCU-HgO  system  in  the  same  packed  column. 
The  values  of  kj-a  obtained  by  Chu  are  18  percent  higher  than  those  obtained 
in  this  investigation.  A  comparison  plot  is  given  on  Plate  VII. 

CONCLUSION 

The  results  of  the  present  investigation  lead  to  the  following  con- 
clusions : 

(1)  This  investigation  fails  to  show  any  different  effect  of  flow 
types,  conntercurrent  and  co-current,  on  the  value  of  liquid-film  co- 
efficient within  the  range  of  study.  The  significantly  higher  liquid- 
film  coefficient  found  by  Collins  (2)  in  wetted-wall  column  for  the 
co-current  flow  did  not  appear  in  the  present  investigation. 

(2)  Discrepancies  between  the  present  results  in  disc  column  and  those 
quoted  in  the  literature  have  been  found.  This  inconsistency  also  exbrta 
among  other  investigator's  work. 

(3)  Though  the  performance  is  quite  different  from  one  disc  column  to 
another,  the  data  from  the  same  column  are  likely  self-consistent,  as 
justified  by  the  agreement  of  the  predicted  correlation  for  Clg-H^O  system 
with  the  experimental  results. 

(4.)  The  empirical  correlation  for  liquid  fi3m  coefficient  in  the  disc 
column  was  found  to  be: 

kL  »  0.0203  p°*745 
This  correlation  is  better  than  both  Hwu's  and  Chu's  correlations  in  view 
of  the  successfulness  in  predicting  CI2-H20  system  data. 
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(5)  A  correlation  of  liquid  film  coefficient  with  liquid  rate  for 
Cl^-H-O  system  in  the  disc  column  is  suggested  as: 

^  -  0.0178  p  °*7^5 
or  the  generalized  equation  can  be  represented  by: 

D 

(6)  The  gas-film  coefficient  for  NH^-HgO  system  in  the  disc  column 
can  be  represented  by: 


V  __0.0326r°-30  (^.33(^-)-0.56  ^ 


T9d 


(7)  The  liquid  film  coefficient  for  the  packed  column  is  given  by: 

kja  -  0.0655  L0'85 
or  by  the  general  form: 

kLa  m   84.0  fcfc-f*  (p0-* 

D 
or  expressed  in  transfer  unit  as: 

H^  =  0.103  L0-15 

(8)  The  satisfactory  interpretation  of  CLj-^0  system  data  justifies 
the  relation  given  by  the  Penetration  Theory,  that  kj,  is  proportional  to 
D0'5  rather  than  D.  The  latter  is  predicted  by  Hatta's  film  theory. 
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APPENDIX 


I.  Principal  Constants  of  Columns 

(a)  Disc  column 

Somber  of  discs  35 

Disc  diameter  1.5  cm.  (0.0492  ft.) 

Disc  thickness  0.48  cm.  (0.0304  ft.) 

Dry  surface  area  of  discs  0.218  sq.  ft. 

Mean  perimeter  for  liquid  flow  0.127  ft. 

Tube  internal  diameter  0.0938  ft. 

(b)  Packed  column 

Size  of  packing  8  mm. 

Height  of  bed  5-l/2  inches 

Tube  internal  diameter  2  inches 

Cross  section  area  0.0218  sq.  ft. 

Volume  of  bed  0. 00954  cu.  ft. 
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II.  Preliminary  Experiments 


Test  for  filtration  effect. 

1  test  run  was  performed  in  the  disc  column  tinder  ccuntercurrent 
flow  condition.  The  recorded  data  are  given  as  follows: 


Barometric  pressure 

731  mm. 

,  Hg 

Gas  rate 

5.42  cu.ft./hr. 

Gas  tempera ture,  inlet 

85©F 

outlet 

85°F 

Column  gauge  pressure 

90  mm. 

water 

Liquid  rate 

14.6  lb./hr. 

Liquid  temperature, 

inlet 

30.5°C 

outlet 

ao.5°c 

Sample  So. 

1 

2 

3 

4 

5 

6 

Sampling  time,  min. 
after  run  started 

8 

10 

12 

K 

16 

18 

Filtration 

+ 

Ml 

+ 

- 

+ 

■• 

Apprcx.  height,  cm., 
nozzle  to  filtrate 

5 

11 

9 

Time  of  filtration, 
minute 

6 

12 

10 

Titration  time,  min. 
after  filtration  started 

6 

17 

10 

CO-  absorbed 

Ib./cu.ft.  0.0431  0.0438  0.0445  0.0408  0.0434  0.0431 

Sampling  flasks  were  rubber-stopped  before  filtration. 

The  various  factors  indicated  above  were  so  combined  that  all  the  effects 

in  a  single  sample  were  additive  and  easy  to  detect.  The  result  favors  the 

statement  that  there  is  no  effect  of  filtration  on  CO  absorbed. 

Test  for  Effect  of  Exposure  Time  on  the  Twmnei  Saspling  Solution. 

Three  samples  were  used  for  this  test.  Each  sample  contained  40ml. 


0.051  N  NaOH  and  20  ml.  0.2N  BaCl  solutions,  and  was  put  in  a  400  ml. 

beaker.  Then  the  following  data  were  obtained: 

Sarnnle  No.  12       3 

Exposure  time,  min.        0       15      30 

CO2  absorbed 

mole/initial  mole  0»    0.865     1.37 

of  BfaaH,  xl04 

•The  effective  concentration  of  NaOH  was  determined 
under  this  assumption . 

Ho  correction  for  COg  so  absorbed  has  been  made  in  evaluating  kr. 

Test  for  Time  for  Reaching  Steady  State  in  The  Packed  Column. 

A  test  run  for  this  purpose  was  performed  under  co-current  flow 

condition.  Observed  data  and  results  are  given  below: 

Barometric  pressure        752.8  mm.  Hg.   Room  Terap.  93°F 

Gas  Rate  5.42  cu.  ft./hr. 

Gas  temperature,   inlet    83°F       outlet  83°F 

Liquid  rate  42.9  Ib./hr. 

Liquid  temperature,  inlet    27.8°C     outlet  27.9°C 

Sample  No.        1        2       3       4        5 

min.  ar  ter 

run  started       4        6       8       10       12 

CO2  absorbed 
Ib./cu.ft.     0.0218     0.0236    0.0233    0.0221     0.0236 

The  result  of  this  run,  though  quite  inconsistent  with  other  runs, 

shows  that  the  data  taken  from  the  4th  min.  deviate  from  the  mean  within 

5  percent.  It  may  also  be  noted  that  the  small  standard  deviation  of  the 

result  of  RE  series  (for  the  packed  column)  leads  to  the  same  conclusion. 


III.  ibeperimental  Data 

Sample  Calculation. 

(a)  Observed  data  for  Run  R00-1A,  co-current  flow,  packed  column, 

Barometric  pressure  744. 8  mm.  Hg  at  75°P 

Gas  rate  1.55  cu.ft./hr. 

Gas  temperature,  inlet  78°F    outlet  7S°F 

Column  gauge  pressure  2.2  cm.  water 

Liquid  rate  13.0  lb./hr. 

Liquid  temperatxire ,  inlet  30.8°C  outlet  23.2°C 

COg  concentration  in  water,  at  inlet        0.0000 

at  outlet      0.0397  Ib./cu.ft. 

(b)  Published  data 

Correction  factor  for  800  mm.  brass  scale  barometer  (13)  -0.130  nm./°C 
Vapor  pressure  of  H20  at  78°F,  from  Keenan's  Thermodynamical  Proper- 
ties of  Steam  p. 28,  »  0.9666  in.  Hg,  or  24  mm.  Hg 
Henry  constant  (13)  at  30.8  -  9.44  x  103  m.  Hg/(lb.  CC^/cu.ft.) 

at  28.2°C  8.88  x  10?  an.  Hg/Clb.CC^/cu.ft.) 

(c)  Calculation 

p.p.  of  C02  in  the  column  ■  718  mm.% 

Ce  te  Fco2  t*   =  °»°76l  Ib./cu.ft.  at  inlet, and 

■  0.0809  lb./cu.ft.at  outlet 
4C=  (C^)^  -  (Cg-S)^  -  0.0397  Ib./cu.ft. 
3a(Ce-C)in  -  *#4<*)U*  -  °'6l5 

^C)l.m.-^C/  3a(Ce-C)ln  -  3a(Ce-C)out  -  0.0568  Ib./cu.ft. 
M  "  (cout  "  cin)  ^  Liquid  rate  /  62.4.  =  8.41  x  1Q~3  lb/hr. 
A  -  0.0218  sq.  ft. 


Substitute  the  above  values  Into  equation  (l),  and  obtain  the  value 
of  M  at  average  temperature  29.5°C,  15.5.  Correct  this  value  to  25°C 
with  equation  (3)  on  page  fl  ,  and  the  resulting  kLa  will  be  14.1 
(ib./hr.-cu.ft. )(lb./cu.ft. ). 
Data   Experimental  data  are  listed  in  the  following  tables: 

Table  3.  Absorption  data  of  CO^O  system  in  the  disc  column  (PC). 

Table  U.     Li-rdd  film  coefficient  for  CO^O  absorption  in  the 
disc  column  (PC). 

Table  5.  Absorption  data  of  CO^O  system  in  the  disc  column  at 
gas  rate  of  5.42  cu.ft./hr.  HW« 

Table  6.  Absorption  data  of  00^0  system  in  the  packed  co'umn  (RE). 

Table  7.  Liquid  film  coefficeint  for  COj-HjO  absorption  in  the 
packed  column  (RE  j • 

In  recording  gas  rate  the  effect  of  variation  of  temperature  and  pressure 
was  ignored.  This  would  introduce  a  maximum  error  of  less  than  3  percent 
(tf.  Catalog  98-A,  Fisher  &  Porter  Company). 

The  differences  between  gas  inlet  and  outlet  temperatures  were  kept 
within  1?F  ,  and  therefore  only  the  inlet  temperature  was  listed  in  the 

tables. 

Since  the  operating  liquid  temperatures  fluctuated  around  25°C,  all 
observed  data  c?  liquid  film  coefficient  were  corrected  to  25°C  rather  than 
20°C,  in  order  to  reduce  effect  of  any  error  associated  with  the  correction 
equation  (3).  However,  for  comparison  with  published  data,  equations  drawn 
from  data  at  25°C  were  corrected  to  20OC  wich  equation  (3). 
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n 

:.  Statistical  Analysis 

PC  series 

.  Data  taken  from  the  last  column  of  Table  4  are  summarized 

In  Table  8a.  The  following  sums  of  squares  are  calculated  from  the 

main  table  and  sub-tables  of  Table  8  a,b,c,  and  d. 

C.F. 

-  (74.65U)2/72 

=  77.40043  ,  where  72  is  the  total  number  of 
observation. 

(a) 

0.7L  :24-—-rl.34242  -C.F.  =  2.49332 

(b) 

.+  2.60962)/2  -  C.F.  -  2.37931 

(c) 

,r  19.33622)/18  -  C.F.  -  0.22781 

(d) 

.r  39.22282)/36  -  C.F.  =  0.19995 

(e) 

f  5.06262)/4  -  C.F.  »  2.04449 

(f) 

t  9.93472)/8  -  C.F.  -  1.93212 

(g) 

+ 37.99212)/36  -  C.F.  -  0.02467 

(h) 

(9.69*22  t 

t  U.9618*)/12  -  C.F.  =  1.92405 

(i) 

t  24.73202)/24  -  C  F.  =  0.02553 

(J) 

+  12.39822)/12  -  C.F.  -  0.05986 

(k) 

f  29.56562)/24  -  C.F.  -  1.89697 

Sum  of  square  for   1  -  d 

,   T  .  e  ,  RxT  =  b-d-e 

for  main  effects: 

F  -  g,       %  -  k,   •  «-  f 

for  interactions: 

FxL  =  h-g-k  ,  FxG  .  j-i-g  ,  IoeG  =  f-k-i 
FxIotG  =  *~f-g-h-i-j-k 

for  observation  within  run:   a-b 

QD  series. 

No  attempt  was  made  to  test  for  significance  of  effects  flow 

type  with 

this  experiment 

.  The  purpose  of  conducting  thi3  experiment  was  to 

furnish  me 

re  data  to  determine  the  relationship  between  k^  and  the  liquid 

rate.  Therefore  only  error  of  observation  within  run  was  calculated.  This 

was  done  by  the  method  as 

described  for  PC  series,  i.e.  to  find  the  corres- 

ponding  a- 

b  term. 

33 


RE  scries.  Computation  similar  to  that  don©  in  PC  series  was  carried  out. 
There  were  two  missing  data,  EOl^lA  and  E31-1&,  due  to  known  errors.  In 
order  to  faci"*tate  computation,  these  data  were  substituted  by  values 
that  would  give  minimum  errors.  The  new  values  were  calculated  with  the 

folowing  formula  (5): 

"      nX.,  +(v-l)X, .  +  X2.  -  X.. 
X||-  (v-2)(n-l) 


X2  |  = 


nX.,  +(v-l)X2.  +  X, .  -  X.. 
(v_?)(n_1} 


where    X  -  value  of  observation 

A 

X..  *:  the  estimated  value  of  the  missing  one  of  replicate  1  and 
treatment  X. 

121  =  the  estimated  value  of  the  missing  one  of  replicate  1  and 
treatment  2. 

*.l  =  the  stm  of  r9Plica1i?!  1  excluding  X^p 

X^#  ■=■   the  sua  of  treatment  1  excluding  Hy. . 

X>#  =  the  grand  total  excluding  the  missing  one. 

a  =  number  of  replicate. 

v  =  number  of  treatment. 
According  to  the  above  formula  and  using  the  first  sample  (a)  of  replicate 
1  as  the  replicate  in  the  formula,  we  can  calculate  the  best  estimate 
value3  for  E0L-1A  and  E31-1A.  These  are  1.2096  and  1.7221  respectively. 

Since  we  have  introduced  two  values  v/lth  pJuimun  error,  the  total 
number  of  observation  for  determining  experimental  error  will  thus  be 
reduced  by  two.  So  total  degree  of  freedom  becomes  96-1-2  ^93. 
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V.Derivation  of  Equations 
(4r)    kL  *  0.0203  p°'7/<-5 
(S)   kLa  =  0.0655  L0,85 

Selection  of  data.   All  data  from  BE  series  were  used  to  derive  equation 
(  5"),  but  data  from  PC  scries  were  examined  carefully,  because  the 
analytical  result  (Table  l)  showed  the  significance  of  variation  due  to 
replicate  in  PC  series.  There  are  two  possibilities;  (a)  the  experiment 
is  lrreproducible  in  the  disc  column,  or  (b)  some  unnoticed  error  has 
Veen  introduced  in  either  of  the  two  replicates. 

Though  the  performance  in  different  disc  colutms  may  give  different 
results  (2.0)  t   it  is  net  plausible  that  the  experiment  is  irreproducible 
in  the  same  column.  The  changable  de-wctting  phenomenon  will  definitely 
increase  the  experimental  error  to  a  considerable  dogree,  tut  it  will 
not  introduce  any  systematical  error.  Therefore,  the  irreprodttdbility 
of  the  experiment  will  not  be  considered. 

The  detail  of  experimental  procedure  was  exactly  the  came  for  both 
replicates.  The  HC1  solution  used  belonged  to  the  same  batch  and  was 
analysed  from  time  to  time  for  check.  The  only  difference  between  the  rep- 
licates was  the  NaOH  solution  used.  In  replicate  1  a  batch  of  HaOH  solution 
of  concentration  0.0465N  was  used,  and  in  replicate  2  another  bateh  of 
0.0510N  was  used.  Two  saiiiples  of  HaOH  solution  vero  det  r-mined  for 
effective  concentration  for  each  batch,  and  one  check  was  made  during  the 
experiment  (For  3D  and  RE  series,  the  concentration  of  CO-  in  inlet  water 
and  the  effective  concentration  of  NaOH  were  checked  every  six  runs.  The 
three  series  were  performed  in  the  order  of  PC,  QD  and  RE.).  Thus  the 
concentration  determination  would  not  be  the  source  of  error.  Since  the 


70 


record  did  not  show  any  perceivable  mistake,  a  comparison  of  data  from 
replicate  1  and  2  of  PC  series  with  those  from  QD  series  was  made. 

The  total  of  164.  k^  data  for  the  disc  column  were  separated  into 
three  groups,  PC-1  (for  rep.  l),  PC-2  and  QD.  Every  observation  was 
compared  with  the  mean  value  of  the  observations  that  belonged  to  the 
same  group  and  had  the  same  liquid  rate  (ignoring  the  flow  type  and 
gas  rate  difference).  If  that  observation  had  a  deviation  more  than  20 
percent  from  the  mean,  it  was  discarded.  The  rejected  observations  were 
P10-2B,  P12-2B,  P20-3B,  P20-2B,  P22-2A,  C02-1B,  C20-U3,  C21-1A,  C10-1A, 
C01-2A,  C21-2A,  Q6-3A  and  B,  D6-3A  and  B,  a  total  number  of  15  observations. 
The  resulting  means  after  these  observations  were  rejected,  were  plotted 
on  Plate  IX.  It  is  obvious  from  the  plot  that  data  from  PC-1  are  con- 
sistent with  those  from  QD,  while  data  from  PC-2  give  higher  kL  values. 
Therefore,  all  data  from  PC-2  were  not  used  in  correlation. 
Derivation.   Least  square  method  (same  as  that  used  in  linear  regression 
determination)  was  used  in  derivation.  The  procedure  can  be  found  in 
Snedecor's  Statistical  Methods  p. 138. 

The  sample  standard  deviation  from  the  resulting  function  (regression) 
can  be  evaluated  by: 

where  n  is  the  number  of  data  used,  b  is  the  slope  (coefficient  of  re- 
gression), y  and  xy  are  deviation  of  I  (dependent  variable)  and  XT 
(independent  variable  times  X)  of  a  single  point  with  the  means.  Since 
Sy.x  is  sxP^ssed  in  the  logarithmic  scale,  we  must  change  it  back  to 
unit  related  to  kL,  e.g.  a  0.02  unit  of  S^x  has  an  anti-logarithm  of 
1.047.  It  means  that  the  standard  deviation  is  4.7  percent. 

Data  beyond  95  percent  confidence  limits  were  rejected.  For  the  disc 
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column  there  were  20,  and  for  the  packed  column  8.  The  remaining  data 
were  used  to  recalculate  the  equation  of  least  deviation.  The  resulting 
equations  were  further  transferred  to  that  for  20°C  by  use  of  equation 

(3). 

The  final  equation  for  the  disc  column,  based  on  PC-1  and  QD,  is: 

^  -  0.0203  V°'U5  to 

and  for  the  packed  column  is: 

kL.  -  0.0655  L0-85  (5) 

and  for  the  disc  column  based  on  PC-2  is: 

T^   -  0.0304  r0,738  (18) 

The  last  equation  given  is  for  comparison  only.  Equation  (4)  is  the 
sole  one  that  is  considered  as  the  result  c?  the  present  experiment  for  the 
disc  column,  and  used  anywhere  for  interpretation  and  correlation.  The 
fact  that  the  two  lines  expressed  by  equations  (4)  and  (13)  are  parallel 
indicated  that  the  experimental  error  associated  with  k^  is  in  simple 
multiple  form  (cf.  p. |I ).  Therefore  use  of  logarithmic  transformat5.on  is 
justified. 

VI.  Derivation  of  Equations 
(8)  and  (9) 

Equations  (8)  and  (9)  are  of  the  general  forms  for  the  disc  column  and 

packed  column  respectively.  They  are  easy  to  obtain  by  comparing  the 

equations  (8)  and  (9)  with  the  generalized  equation  (6)  and  (7)  together 

with  the  following  data: 

JL-    2.42  lb.mass/hr.ft.   (13) 

^=    62.4  lb.  mass/cu.ft.     (13) 
X>  =    7.0  x  10-5  ft2/hr.,       (17) 


VII.  Equations  for  Clg-HgO 
System  in  the  Disc  Column 

Eauation  nredicted  from  the  result  of  this  investigation.  From  the 

generalized  equation  ( L  ) 

,  kL  is  proportional  to  D   .  The  equation  for 

C02-H  n  system  is: 

kL=  0.0203  P0,745 

and  the  diffusivity  of  Clg  in  water  is  5.4  x  10~5  ft2/hr.  %&  ).  There- 

fore, the  equation  for  CLj-B^0  wm  °«* 

k.=  0.0178  p0'14-5                                             (12) 

Equation  predicted  by  use 

of  Ghu's  eauation.   Chu's  111  eauation  for 

CCL-EjO  system  in  the  same  disc  column  is: 

kL=  0.03R3r°-S44 

and  the  resulting  equation  for  CLj-H-O  system  is: 

kL=  COS^Sf0,644      .            (19) 

Some  particular  points  or 

i  lines  predicted  by  the  eauations  (  )  and  (  n), 

For  eouation  (is ): 
For  equation  (i2)t 

r          r          $ 

150      25.1     0.848 
200      30.3     1.03 
300      39.3     1.33 
400      47.5     1.61 

r           p'^     ,,, 

150      42.0     0.748 
200      50.5     0.897 
300      70.0     1.240 
400      87.0     1.545 

.  J 


VIII.  Gas-Filn  Coefficient  of 
Ammonia-water  System  in  the  Disc  Column 

Equations* 

(a)  Predicted  equation  for  liquid-film  coefficient* 

The  diffusivity  of  NH^  in  water  is  given  as  7.9  x  10"*5  ft.2/hr. 
(19).  The  equation  for  CO2-H2O  system  in  the  disc  column  was  suggested 
by  the  present  paper  as: 

i^  -  0.0203  rl0#745 
and  so  the  resulting  equation  for  HHy^O  system  will  be  (ef .  Appendix  VII ) : 

1^=  0.0216  P°'U5  (20) 

(b)  Relation  between  overall  and  film  coefficients  is: 

~I_  =    -I*.  *  JL.        (21) 

*g    h  \ 

as  given  in  any  absorption  book.  The  fienry  constant  H  for  amnonia  in 
water  at  20°C  is  taken  as  0.013  atm./(lb.  of  MHycu.lt.  of  H20),  which  was 
evaluated  by  Hwu  (9)  from  Kowalke's  equation. 

(c)  Jcq  vs.  relative  velocity: 

Though  Hwu's  experimental  result  showed  that  the  relation  between 

kg  and  relative  velocity  in  the  disc  column  could  be  represented  by 

0  6? 
kg^v  *  ,  the  following  conventional  relation  was  used: 

V-v0-6?  (22) 

which  has  been  verified  by  Stephens  and  Morris,  and  also  by  Taylor  and 
Roberts  in  four  different  disc  columns  (19)  (20). 


Experimental  Data. 

The  RH  absorption  experiment  in  the  disc  coltaan  was  performed  by 
Hwu  (9).  Data  at  six  d5.fferent  liquid  rates  with  nearly  equal  intarval 
were  taken  for  the  present  correlation.  These  data  and  the  corresponding 
k-  values  calculated  with  equations  (20)  (2l)  and  (22)  are  listed  in  the 
following  table: 

Table  11.  Absorption  data  of  NH  -H,-0  system 


Run  :  Liquid 
No.  :   rate   : 

:  Relative 
velocity 
:  ft. /sec. 

i   ^ 

: 

at  20°C 
lb. 

• 
• 

*G 

:lb./hr.ft. 

:   hr 

.sq.ft.atm. 

:   At 

V 

5.84  ft./s« 

2      220 

A.  82 

9.55 

10.6 

12.1 

3      360 

5.34 

12.2 

13.4 

L4.2 

5      185 

4.69 

8.91 

10.0 

11.6 

7      255 

4.97 

10.57 

11.8 

13.2 

10      300 

5.33 

11.6 

12.9 

13.7 

11     395 

5.65 

13.02 

14.3 

14.6 

The  first  four  columns  were  taken  from  Table  3  of  Hwu* s  thesis  (9). 
The  last  two  columns  were  calculated  with  equations  (20)(2l)  and  (22). 

With  the  data  given  above  we  can  derive  an  equation  of  kg  vb, 

through  the  same  procedure  as  described  in  Appendix  V.  The  resulting 

Squation  is: 

kG  =  3.06P°-30         (23) 

In  order  to  compare  with  other  investigator's  results,  we  transform 

equation  (23)  to  satisfy  the  condition  of  relative  velocity  equal  to 

8.4  ft. /sec.  Thus  the  fanal  equation  becomes: 

k^  -3.90P0'30 

At  two  particular  points,  p^  155  Ib./hr.ft.    kj,  ^17.7  lb./(hr,  ft?atm.) 

P  =400  lb./hr.ft.    kg  -  23.6  lb./(hr.ft2.atm.) 
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ABSTRACT 

The  purpose  of  this  investigation  was  to  study  the  effect  of 

flow  type  on  liquid  film  coefficient  in  gas  absorption.  Experiments 

under  both  countercurrent  flow  and  co-current  flow  were  carried  out, 

and  their  results  were  compared.  Experimental  data  of  some  previous 

workers  were  also  reeorrelated. 

CO  -H„0  system  was  chosen  for  this  study.  Apparatus  used  were 
2  2 

1-2/2  inch  ^sc  column  with  35  pieces  of  1.5  mm  diameter  X  0.4.8  th»k 
ceramic  discs,  and  a  2  inch  packed  column  with  3  mm  packings.  The 
experiments  were  designed  end  their  results  were  analyzed  statistically. 
The  liquid  rate  covered  a  range  from  10.0  lb./hr.  to  52.6  lb./hr.,  and 
gas  rate  covered  a  range  from  1.55  cu.ft./hr.  to  10.6  cu.ft./hr. 

From  the  results  of  this  investigation  the  following  conclusions 
were  reached: 

(1)  Flow  type  has  no  effect  on  magnitude  of  liquid  film  coefficient 
within  the  range  of  .study. 

(2)  Results  obtained  from  disc  column  will  generally  vary  from  column 
to  column. 

(3)  Data  of  COg-ikj^  system  for  the  columns  used  in  this  investigation 
can  be  correlated  by  the  following  equations: 

for  the  disc  column  l<L  =  0.0203  p0,7^5  at  20°C, 

0  3*5 
for  the  packed  column        kLa  =0.0655  I»  *   at  20°C. 

(4.)  A  general  equation  for  liquid  film  coefficient  in  the  disc  column 

can  be  obtained  from  the  present  results.  It  is  expressed  asx 


This  equation  was  verified  with  experimental  data  of  C12-BL0  system 
obtained  by  the  previous  investigators. 

(5)  Gas  film  coefficient  in  the  disc  column  can  be  predicted  with 
considerable  accuracy  with  the  following  equation: 


a?   ^  0.0326  p°.30  x  (^)-0.33  (-^)-0.56   _p_ 
obtained  f n 
another  investigator. 


which  was  obtained  from  the  results  of  this  investigation  and  that  of 


